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Abstract 

Approximately 438 pb _1 of e + e~ data from the OPAL detector, taken with the 
LEP collider running at centre-of-mass energies of 192-209 GeV, are analyzed to 
search for evidence of chargino pair production, e + e~ — ► XiXi> or neutralino asso- 
ciated production, e + e~ — > X2Xi- Limits are set at the 95% confidence level on the 
product of the cross-section for the process e + e~ — ► XiXi an d its branching ratios 
to topologies containing jets and missing energy, or jets with a lepton and missing 
energy, and on the product of the cross-section for e + e~ — > x^,X\ an d its branching 
ratio to jets. R-parity conservation is assumed throughout this paper. When these 
results are interpreted in the context of the Constrained Minimal Supersymmetric 
Standard Model, limits are also set on the masses of the xf,Xi an d X2> an d regions 
of the parameter space of the model are ruled out. Nearly model- independent limits 
are also set at the 95% confidence level on cr(e + e~ — > XiXi ) with the assumption 
that each chargino decays via a W boson, and on a(e + e~ — > X2X1) with the X2 
assumed to decay via a Z°. 
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1 Introduction and Theory 



Supersymmetric (SUSY) extensions of the Standard Model postulate the existence of 
fermionic partners for all Standard Model bosons, scalar bosonic partners for all Standard 
Model fermions, and at least two complex Higgs doublet fields, each containing a charged 
and a neutral component [1]. The partners of the U(l)y weak hypercharge field B° 
and the neutral component of the SU(2) weak isospin field W° (gauginos), and of the 
two neutral components of the Higgs fields (higgsinos), mix to form the four neutralinos 
X° 2,3,4; where the index increases with the neutralino mass. The higgsino partners of the 
charged components of the Higgs fields and gaugino partners of the charged components 
of the weak isospin field mix to form the charginos xf 2 > where the index 1 indicates the 
lighter chargino. 

SUSY particles couple to the same Standard Model particles as their Standard Model 
partners. The lightest charginos may thus be pair-produced at LEP in s-channel 7 or 
Z° exchange. Their gaugino components can also be pair-produced in t-channel scalar 
electron neutrino (u c ) exchange if the v c is sufficiently light. The chargino pair production 
cross-section is expected to be large [2]. It is typically a few picobarns (pb) and depends on 
the mass and mixing parameters of the SUSY model; however, the interference between 
the s-channel and t-channel processes is destructive, so if there is a light v c the cross- 
section can be small. If R-parity conservation is assumed, as it is throughout this paper, 
the lightest SUSY particle (LSP), which is expected to be either the x± or a ^> is stable and 
undetectable. If the LSP is the Xi, neutralinos are most likely to be detected in associated 
production of X2X15 either through s-channel 7 or Z° exchange or by t-channel e exchange 
if there is a light scalar electron. Neutralino cross-sections are comparable to those for 
chargino pair production only if the higgsino component is predominant; however, if 
charginos are too massive to be pair-produced, neutralino associated production may be 
the only SUSY process with a visible signature at LEP [3]. In such cases, the neutralino 
search can be used to rule out regions of SUSY parameter space which are not accessible 
to the chargino search. A list of the chargino and neutralino production and decay modes 
considered in this paper, along with a few final states which were not considered, is given 
in Table 1. 

In R-parity-conserving SUSY models where the Xi is the LSP, the xt decays to a 
Xi and a W ± : which is off-shell if the mass difference AM± between the xt an d the Xi 
is less than the mass. The W^ 1 in turn decays either to qq or lis, the latter being 
increasingly favoured if AM± is very small. If there is a light u, an additional decay 
channel xt ~^ vft 1 opens up, and will predominate if the v is the LSP. If the v is not 
the LSP, it will decay to x\ v - Neither of these decay channels leads to a new final-state 
topology; they simply increase the fraction of leptonic final states. The existence of light 
scalar leptons would lead to similar leptonic final states. The event topologies for chargino 
pair production contain missing energy from undetected x? or v and neutrinos, and either 
four jets, two jets and a lepton, or two leptons. The analysis of the fully leptonic topology 
has been published separately [4]. The kinematics of the event depend mainly on AM±, 
since that determines the energy remaining for the Standard Model particles in the final 
state. There is also some dependence on the mass of the xt ■> which determines the boost 
of the final state particles. 

Similarly, a X2 i s expected to decay to a Xi an d a virtual Z°, which will decay to 
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Table 1: Principal production and decay modes for the lightest charginos and neutralinos in e + e~ colli- 
sions. All detectable final states consist of missing energy and jets or charged leptons or both. The final 
states above the horizontal line in each category are considered in this paper. The chargino three-body 
leptonic decays are analyzed in [4], and the efficiencies, background estimates, candidates and likelihood 
distributions from that analysis are used in this paper to set cross-section limits on chargino production. 
Chargino two-body decays and neutralino invisible and leptonic decays listed below the lines are not 
considered here when setting limits, although the chargino two-body leptonic decays with the stable v in 
the final state are analyzed in [4]. 



Chargino pair production 


e + e-^xrxr^(W + *^)(W-*^)^ + ^'-^Xi 

->qq'Xi&4Xi 
-qq'XiWXi 


^(u£ + )(u'£'-) ^stable (z> LSP) 


Neutralino associated production 




^i~x\x\ 

-^XiX? 



qq, £ + £~ or vv. In the presence of light scalar leptons or scalar neutrinos, additional 
decays become possible which increase the £ + £~ or the invisible vv branching fractions. 
Detectable event topologies consist of missing energy from the two undetected Xi particles 
and a pair of jets from qq or a pair of leptons. The kinematics depend on AM , the mass 
difference between the x° an d Xii i n the same way as for charginos and on the mass of 
the X2- If AMo is small, the experimental signature becomes monojet-like because of the 
small angular separation between the qq jets. 

The analysis for chargino and neutralino decays into hadronic and semileptonic final 
states described in this paper is a new one, different from the one used to analyze the 
OPAL data at centre-of-mass energies below 192 GeV [5-7]. The results for charginos 
decaying to leptonic final states are taken from the analysis described in [4]. 

2 Detector and Data Samples 

A complete description of the OPAL detector can be found in [9] and only a brief overview 
is given here. The central detector consisted of a system of tracking chambers covering 
the angular 1 region |cos#| < 0.97, inside a 0.435 T uniform magnetic field parallel to 
the beam axis. It was composed of a large-volume jet chamber surrounded by a set of z 
chambers measuring the track coordinates along the beam direction and containing a high 
precision drift chamber and a silicon microstrip vertex detector [10]. Scintillation counters 
surrounded the solenoid and scintillating tiles [11] were also present in the endcaps just 
outside the pressure bell for the central tracking chambers; these provided time of flight 

1 The OPAL coordinate system is defined so that the z axis is in the direction of the electron beam, 
the x axis points toward the centre of the LEP ring, and 9 and <j> are the polar and azimuthal angles, 
defined relative to the +z and +x axes respectively. 
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information. A lead-glass electromagnetic calorimeter, the barrel section of which was lo- 
cated outside the magnet coil, covered the full azimuthal range with excellent hermeticity 
in the polar angle range of | cos#| < 0.82 for the barrel region and 0.81 < | cos#| < 0.984 
for the endcap region. The magnet return yoke was instrumented for hadron calorimetry 
and consisted of barrel and endcap sections along with pole tip detectors that together 
covered the region | cos#| < 0.99. Four layers of muon chambers covered the outside of 
the hadron calorimeter. In the region |cos#| > 0.99, electromagnetic calorimeters close 
to the beam axis completed the geometrical acceptance down to about 33 mrad. These 
included the forward detectors, which were lead-scintillator sandwich calorimeters and, 
at smaller angles, silicon tungsten calorimeters [12] located on both sides of the interac- 
tion point and used to monitor the luminosity, and the forward scintillating tile counters. 
The gap between the endcap electromagnetic calorimeter and the forward detector was 
instrumented with an additional lead-scintillator electromagnetic calorimeter, called the 
"gamma catcher" . 

The data taken with the OPAL detector in 1999 and 2000 at y/s = 192 - 209 GeV 
were studied for the hadronic and semileptonic channel analyses described in this paper 
(see Table 2 for details). Data were used for this search only if the central jet chamber, 
barrel and endcap electromagnetic calorimeters, hadron calorimeters and forward detec- 
tors, including the silicon tungsten calorimeters, were all functioning. For convenience, 
Table 2 also shows which data were studied in the leptonic channel analysis [4] , which is 
combined with the other two channels to produce limits on chargino pair production. 

Table 2: Data samples used in the analyses described in this paper (qq for x°Xi associated production 
and qq^, qqqq for chargino pair production) as well as in the lepton pair analysis (££) [4]. Energy bins 
are defined so that v^min < ^ V^max- The total luminosity studied is 437.6 pb -1 in the semileptonic 
and hadronic channels and 440.4 pb -1 in the leptonic channel. 





qq£, qqqq 




t 


t 




V^min 


V^max 


V^avg 

(GeV) 




V^min — V^max 


V^avg 

(GeV) 


!C 


(GeV) 




(Pb- 1 ) 


(GeV) 


(Pb" 1 ) 


190.0- 


194.0 


191.6 


29.14 


190.0 - 194.0 


191.6 


29.3 


194.0 - 


198.0 


195.5 


73.96 


194.0 - 198.0 


195.5 


76.4 


198.0 - 


201.0 


199.5 


75.40 


198.0 - 201.0 


199.5 


76.6 


201.0 - 


202.5 


201.7 


38.27 


201.0 - 204.0 


202.0 


45.5 


202.5 - 


204.5 


203.7 


10.08 








204.5 - 


205.5 


205.1 


71.94 


204.0 - 206.0 


205.1 


79.0 


205.5 - 


206.5 


206.3 


65.56 








206.5 - 


207.5 


206.6 


64.98 


206.0 - 208.0 


206.5 


124.6 


207.5 - 


210.0 


208.0 


8.25 


208.0 - 210.0 


207.9 


9.0 



Chargino and neutralino signal events are generated with the DFGT generator [13], 
which includes spin correlations and allows for a proper treatment of the and Z° boson 
width effects in the chargino and neutralino decays. The generator includes initial state 
radiation (ISR) and uses JETSET 7.4 [14] for the hadronization of the quark-anti-quark 
system in the hadronic decays of charginos and neutralinos. Signal events are generated 
at centre-of-mass energies of 196, 202 and 207 GeV. For chargino pair production, 5000 
events are generated at each point on a grid with 5 GeV intervals in chargino mass starting 
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at 75 GeV, with the last point at 0.5 GeV below half the beam energy. For example, for 
= 196 GeV, the generated chargino masses are 75, 80, 85, 90, 95 and 97.5 GeV. For 
each chargino mass, 22 equally spaced lightest neutralino masses are generated, from zero 
to 3 GeV below the chargino mass. For X2X1 production, 5000 events are generated at 
each point on a grid with 2, 5, 10 or 20 GeV spacing in AM from 3 GeV to M(%°), with 
closer spacing for smaller AM , and 5, 10 or 20 GeV spacing in M(xV) + M(x1) from 
100 GeV up to 1 GeV below the centre-of-mass energy, with closer spacing nearer to the 
centre-of-mass energy. 

The main sources of background (see Table 3) are events with genuine missing energy 
either from neutrinos, as in the case of four-fermion final states including a W — > Iv or 
Z° — > vv decay or any events with r decays, or from particles escaping down the beampipe 
as in events with ISR and "two-photon" interactions, where the interaction is between 
initial state photons radiated by the e + and e~ (e + e~ — > e + e~77 — > e + e~ //) rather than 
directly between the electron and positron. In referring to two-photon events, "untagged" 
events are those where the photons are both real and the electron and positron are lost 
in the beampipe, "singly-tagged" means that one photon is sufficiently off-shell to kick 
an electron into the detector, and in "doubly tagged" events, both photons are virtual 
and are seen in the detector. The generated luminosities are generally several hundred 
times the data luminosity, although for the two-photon processes they were only a few 
times to a few tens of times the data luminosity, and the available Bhabha luminosity at 
205-207 GeV was only a few times the data luminosity. 

Signal and background events are processed through a detailed simulation of the OPAL 
detector [27] and analyzed in the same way as the OPAL data. 

3 Analysis 

The main goal of this analysis is to have a selection which is optimized at every kinemat- 
ically accessible value of (M-±,AM±) or (M^>, AM ), henceforth shortened to (M, AM) 
when referring generically to the chargino and neutralino analyses. This ensures that 
there are no "efficiency valleys" , values of (M, AM) where a signal might be missed be- 
cause the analysis was optimized for other regions. In order to make this task feasible, 
the same variables must be used for all (M, AM) so that the optimization can be done 
systematically. Since many of the most useful selection variables are properties of lep- 
tons or jets, a single set of variables cannot describe all of the possible signal topologies. 
The chargino pair production (xtxi) analysis is therefore split into three channels: fully 
hadronic, semileptonic (one lepton and some jets) and fully leptonic. 

The analysis of fully leptonic final states has been published separately [4]. This 
paper describes the analyses of the hadronic and semileptonic final states, which arise 
from chargino decays via a W. If all charginos decay to a lightest neutralino and a W, 
which decays to all leptons and to the four lightest quarks when kinematically accessible, 
then the fractions of events falling into each channel will be the same as for W-pair decays: 
46% qqqq, 44% qqlv, 10% l+t'vv. The branching ratios are very sensitive to AM when 
AM is below about 4 GeV. 

The neutralino associated production (X2X1) search is performed only in the fully 
hadronic channel, as this final state is expected to account for 70% of the signal rate if 
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Table 3: Standard Model background cross-sections at -Js = 208 GeV. With the exception of the PHOJET 
and Vermaseren samples, all background samples are generated at all of the centre-of-mass energies in 
Table 2. For processes where a sample was not generated at 208 GeV, the cross-section shown here was 
re-scaled from the 207 GeV sample. The last column shows the approximate cross-section surviving the 
preselection cuts described in Section 3.2. 



Background Process 


Monte Carlo 
Generator 


Cross-se 
Generator 


ction (pb) 
Preselection 


z/z/7 


NUNUGPV [15] 


7.90 




e + e~ (e + and e~ scatter > 12.5°) 


BHWIDE [16] 


493 


15 


e + e~ (near beam axis) 


TEEGG [17] 


637.6 


20 


T + T - 

qq(7) , 


> 


KK2f [18] 


6.815 
6.62 
78.44 


0.2 
0.2 

3 


77(7...; 




RADCOR [19] 

L J 


19.818 




mi 1 
qqqq J 


(four-fermion final states 
without e + e~) 


KoralW a [201 


3.416 
9.426 
8.966 


10 


e + e~e + e~ 
e + e~ pi + 
e + e~r + T~ 
e + e~qq 


(excluding 
> multiperipheral 
two-photon process) 


grc4f [21] 


16.749 
11.575 
1.726 
37.788 


2 


e + e~e + e~ 


two-photon) 


Vermaseren 6 [22] 


876 


0.3 


e + e~ 

e + e~T + r~ 


> (two-photon) 


BDK C [23] 


619.3 
469.0 


0.2 
0.9 


e + e qq (untagged two-photon) 


PHOJET d [24] 
(PYTHIA [25] for 
hadronic channel 
analyses) 


11166 


6 


e + e~qq (single-tagged two-photon) 


HERWIG [26] 


327.8 


0.7 


, (double-tagged 
e e qq v 

two-photon) 


PHOJET 


3.34 





a Four-fermion events without an c + c~ pair in the final state are simulated with KoralW using matrix 
elements from grc4f to ensure a correct treatment of interference among the diagrams, as the treatment 
of ISR in KoralW is superior to that in grc4f. 

6 No Vermaseren samples were generated at 192 GeV, 204 GeV, 205 GeV or 207 GeV. 

Tncludes effects of additional ISR compared with Vermaseren. 

d No PHOJET samples were generated at 192 GeV, 204 GeV, 205 GeV, 207 GeV or 208 GeV. 
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Z branching ratios are assumed. This assumption is valid for AM above about 10 GeV; 
for smaller AM, the leptonic and invisible branching fractions become relatively more 
important as mass thresholds are crossed. 

In order to combine the results of the three chargino channels, they are designed to 
select non-overlapping subsets of the data. There is, by construction, no overlap between 
the hadronic and semileptonic channels. There remains a small overlap between the 
leptonic channel and the others. The overlap has been checked for candidates in the 
signal Monte Carlo passing stringent signal-selection requirements, which dominate the 
limit-setting, and found to be negligible for the hadronic channel and less than 1% for the 
semileptonic channel. 

3.1 Construction of the Likelihood Function 

The description that follows is for the hadronic and semileptonic channels only, although 
the analysis of the leptonic channel is very similar [4]. 

After the common preselection described in Section 3.2, a set of additional cuts, de- 
scribed in Sections 3.3 - 3.5, is applied in each channel. After these cuts, a set of dis- 
criminant variables is computed. The distributions of these variables are formed for each 
of the Standard Model background Monte Carlo samples used, as well as for the signal 
Monte Carlo samples for all the generated (M, AM) points. Much of the background can 
be rejected by rejecting events which do not fall approximately within the same range in 
each discriminant variable as the signal. The ranges of the signal sample variable distri- 
butions are checked at each generated value of (M, AM), as in general they vary with 
(M, AM). A set of fairly loose cuts which reject values of discriminant variables outside 
the signal range is stored for each (M, AM). 

The variables which have the greatest discriminating power against the backgrounds in 
trial likelihoods evaluated at several different (M, AM) values are then used to construct 
a likelihood discriminant. Reference histograms are made for each variable x« for the 
signal and the sum of all major backgrounds after the cuts for each (M, AM) point. All 
histograms are normalized to unity, after which the histogram contents are denoted by 
Ps(xi) and Pb{xi) for signal and background reference histograms. These histograms 
are given as inputs to a calculation [28] which constructs the signal and background 
likelihoods L s HiPs( x i) an d L B ps Ili-PsO^); the likelihoods are not exact products 
of the probabilities because correlations between the variables are projected out [28]. A 
likelihood discriminant L R = L S /(L S + L B ) is then constructed at each (M, AM) for 
every event in the data sample as well as for the signal and background Monte Carlo 
events. 

The same likelihood function is used for the data collected at all centre-of-mass en- 
ergies. Variables are scaled to the centre-of-mass energy wherever possible in order to 
minimize energy dependences, except in the case of four-fermion backgrounds where fea- 
tures of some variable distributions depend on absolute thresholds. The signal samples 
used to construct the reference histograms are generated at = 207 GeV, with a cross- 
check set made at 196 GeV to test the validity of the assumption of energy independence. 

Likelihood reference histograms are made for every (M, AM) point at which signal 
was generated. The generated 5 GeV spacing in M± is adequate, since the kinematics do 
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not change rapidly with M±; however, for AM, a spacing of 2 GeV is required. Since the 
chargino signal samples are generated with AM spacings ranging from about 3.4 GeV for 
M = 75 GeV to 4.8 GeV for M = 103.5 GeV, this involves interpolating signal reference 
histograms for the likelihood. This is done using the linear interpolation technique de- 
scribed in [29,30] and also used in [4]. Similarly, in the case of the neutralino analysis, 
reference histograms are made for every (M, AM) point at which signal was generated, 
and are interpolated to obtain a grid with a 5 GeV spacing in M^o and a 2 GeV spacing 
in AM, which is again found to be adequate. 

3.2 Common Preselection (Charginos and Neutralinos) 

The common preselection used for the hadronic and semileptonic channels of the chargino 
and neutralino searches is designed to select well-measured events with missing transverse 
energy and to eliminate events which are due to cosmic rays passing through the detector 
or events where energy is missing along the axis of the beam due to ISR or escaping initial 
state electrons in two-photon interactions. Events with jets or leptons detected very close 
to the beam axis are rejected since they tend to be poorly measured. This cut also rejects 
a fraction of the events from LEP machine backgrounds, such as those from beam-gas 
or beam-wall interactions. Energy and momentum information from tracks in the inner 
detector and energy deposit clusters in the calorimeters are combined using an algorithm 
which matches tracks and clusters to avoid double-counting [7] . This algorithm is used in 
the construction of all event quantities such as the visible energy, E V - 1S , which is the sum 
of the energies associated with all tracks and clusters in the event, corrected for double 
counting, the missing transverse momentum, p™ lss (which is the component of momentum 
in the x — y plane required to balance the vector sum of the momenta of all the charged 
tracks), the invariant mass of the system comprising all the detected momentum and 
energy, M vis , and the hadronic mass, Mh a d, which is similar to M vis but excludes tracks 
and clusters associated with lepton candidates. The common preselection consists of the 
following requirements: 

• less than 2 GeV may be deposited in any of the forward detectors, including the 
silicon tungsten calorimeters, and less than 5 GeV may be deposited in any of the 
"gamma catchers" ; 

• events are vetoed if they contain activity in the endcap and forward detectors consis- 
tent with forward-going muons, or if the energy deposits in the forward scintillating 
tile counters are large enough to be consistent with a photon, or if the fraction of 
visible energy in the forward region |cos#| > 0.9 exceeds 20% of y/s, or if they 
exhibit evidence of instrumental noise in the jet chamber; 

• there must be enough tracks, excluding any lepton candidates, to group the tracks 
and clusters in the event into two jets using the Durham algorithm [31] and still 
have at least one charged track in each jet; 

• at least 20% of the charged tracks must pass the "good track" criteria [32]; 

• in events with fewer than 5 "good" charged tracks, the sum of the charges of good 
tracks in each jet must be between —1 and +1 and the sum of charges of all good 
tracks in the event must be zero; 
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• probable cosmic ray interactions are vetoed if any of the following conditions are 
satisfied: 

— time-of-flight information is recorded, and differs by more than 10 ns from that 
expected for collision events; 

— time-of-flight is not recorded, there is at least one hit in the muon chambers, 
only two good tracks are present, and their point of nearest approach to the 
interaction point is not within 10 cm; this rejects cosmic rays going through 
the tracking chambers; 

— there is at least one hit in the muon chambers, there are two or fewer good 
tracks, and ten or more clusters in the electromagnetic calorimeters, and clus- 
ters in the hadron calorimeters; this rejects cosmic rays going through the 
calorimeters; 

— there is no hit in the time-of-flight system, and both jets (when the event 
is forced into two jets) have polar angles with respect to the beam axis of 
Icos^jctl < 0.8; this rejects out-of-time events where a particle goes through 
the barrel of the detector; 

• unmodelled low-energy two-photon backgrounds and some poorly modelled forward 
Bhabha events are removed by requiring p™ lss in excess of 4 GeV, M v ; s in excess 
of 5 GeV and acoplanarity (defined for an event forced into two jets as 180° minus 
the angle between the two jet thrust vectors projected into the x — y plane) greater 
than 5°; 

• there must be significant missing energy: the magnitude of the vector sum of charged 
track momenta and the scalar sum of the energy of all electromagnetic calorimeter 
clusters must each be less than 65% of y/s. 

After this preselection, the accepted cross-sections are reduced to the levels shown in the 
final column of Table 3. Signal efficiencies after the preselection vary from about 8% for 
AM = 3 GeV to about 60-70% or more for AM > 10 GeV. 

3.3 Chargino Pair Production Hadronic Selection 

The event topology for the hadronic channel of the XiXi search consists of jets and 
missing energy, the amount depending on AM. There are four primary quarks in this 
final state, but the jets may be highly boosted and difficult to separate, or additional jets 
may be present from gluon radiation. The principal backgrounds are therefore events with 
hadronic activity and missing energy: two-photon e + e~qq, qq(7) in which one or more 
ISR photons are lost down the beampipe or the jets are mis-measured, W + W~ — > qqqq 
where the jets are mismeasured, W + W~ — ► qqiu where the lepton is not reconstructed 
and Wez/ — > qqez/ where the electron goes down the beampipe. 

To avoid overlap with the semileptonic channel, events are vetoed if they contain a 
lepton candidate which would be accepted in that channel. The overlap between the 
hadronic channel and the fully leptonic channel was checked and is negligible. 

To reduce some of the main backgrounds, the following additional preselection cuts 
are made: 
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events are required to have at least five charged tracks; 



• events are removed if the missing momentum vector points to the forward region, 
I c °s#miss| > 0.9; they are probably due to ISR or two-photon interactions or may 
be missing part of a jet down the beampipe; 

• the missing transverse momentum, j9™ lss , is required to be in excess of 6 GeV, 
whether or not information from the hadronic calorimeter is used in the calcula- 
tion; 

• the Durham jet finder must find between three and five jets with a jet resolution 
parameter of y cut = 0.005; 

• when the event is forced into two jets by adjusting the resolution parameter y cut : 

— each must satisfy the requirement | cos6>j Ct | < 0.9 on its polar angle with respect 
to the beam axis; 

— to reduce qq events, the acoplanarity of the jets, acO p, is required to exceed 
10°; 

— the energy of the more energetic jet, £j e ti, must be less than 100 GeV; 

• the invariant visible mass, M vis , must be less than 180 GeV to reduce the four- 
fermion background; 

• if there is a "loose" lepton candidate, an isolated track which does not satisfy the 
requirements of the semileptonic channel (see Section 3.4), its energy, _E loosc _^, must 
be less than 50 GeV. 

As described in Section 3, cuts are made on the likelihood variables (see below) at each 
(M, AM) in order to eliminate events obviously incompatible with signal. After these 
mass- dependent cuts, efficiencies for hadronic signal events vary from about 1% for AM = 
3 GeV to 40-60% for AM > 5 GeV. Two-photon backgrounds are eliminated for AM 
greater than about 10 GeV and reduced to less than 0.01 pb for smaller AM. Four- 
fermion backgrounds vary from a negligible value for AM < 10 GeV to about 0.7 pb in 
the region dominated by W + W~ background with M w 80 GeV and AM M, while 
qq(7) backgrounds vary from a negligible value for AM < 10 GeV to about 0.005-0.05 pb 
for larger AM. Other backgrounds are negligible. 

The likelihood discriminant is then constructed using the following variables: 

• | cos^ m i SS |, the polar angle of the missing momentum vector, discriminates against 
two-photon events where one or both initial state electrons are lost down the beampipe 
and events with ISR; 

• p™ ss also rejects two-photon and ISR events; 

• M vis reduces four-fermion backgrounds; 

• 1/45; 1/34; 2/23, the values of the Durham jet resolution parameter for which the event 
passes from five jets to four, four to three, and three to two, discriminate against 
two jet events from W + W~ — > qqiu with a poorly reconstructed lepton and qqy; 
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• Egr, the energy of the lowest multiplicity jet (where the number of jets is determined 
by the default Durham resolution parameter), has some discrimination against back- 
grounds containing r leptons; 

• when the event is forced into two jets, aC o P , the acoplanarity angle between the jets, 
discriminates against Standard Model two-fermion events and photon conversions; 

• ^jeti,2, the energies of the two jets, reject four-fermion backgrounds such as W + W~ 
and Wez/; 

• cos^j Ct i 5 2, the cosines of the polar angles of the two jets with respect to the beam 
axis, reject two-photon and ISR events; 

• -^Vioose-£, the number of lepton candidates passing loose requirements, less stringent 
than those of the semileptonic channel, helps to reject W + W~ — ■> qqtv events; 

• E\ oose _£, the energy of the most energetic such loose lepton candidate, also helps to 
reject W + W~ — > qqtv events. 

Several of these variables are substantially correlated, but the effect of these correlations 
is minimized by the projection method used to construct the likelihood, as stated in Sec- 
tion 3.1. Distributions of these variables are shown in Figure 1 after the preselection cuts 
have been applied. Likelihood distributions are shown for the signal, expected background 
and observed data for two pairings of (M, AM) in Figure 2. 

3.4 Chargino Pair Production Semileptonic Selection 

The event topology for the semileptonic channel consists of two jets, a lepton and an 
amount of missing energy which is dependent on AM. The principal backgrounds are 
W + W~ — > qqlv, two-photon processes leading to e + e~qq and e + e"r + r" final states, 
t + t~ with at least one of the taus decaying hadronically, and qq(7) events with a lepton, 
usually fake, in one of the jets. 

Overlap with the hadronic channel is avoided by requiring the presence in the event 
of an isolated lepton candidate with the same isolation criteria as in searches at lower 
energies [7]. Overlap with the leptonic channel is almost eliminated, as for the hadronic 
channel, by vetoing events which pass the selection for W + W~ — » i + vl~v events of the 
W + W~ analysis. 

After applying the general preselection and the veto cuts just described, the following 
additional cuts are made. To remove poorly reconstructed events, it is required that the 
event visible mass, M vis , and hadronic mass, M had , as well as the invariant mass of the 
system consisting of the lepton candidate and the missing momentum vector, M^ m i SS , all 
be less than or equal to the centre-of-mass energy. Cuts are made at each (M, AM), 
as described in Section 3, to limit the accepted variable ranges to roughly the range of 
the signal distribution. These cuts are applied to all of the variables used for likelihood 
distributions, which helps to ensure that the signal reference histograms used in the 
likelihood have data in all bins and a smooth distribution. The same type of cut is 
applied to some additional variables which are not used in the likelihood because they 
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Figure 1: Comparison of data and Monte Carlo distributions of selected likelihood variables for the 
hadronic chargino analysis, shown at preselection level for y/s = 206 GeV. The shaded histograms show 
the various expected backgrounds (the darkest shades are two-fermion processes, predominantly qq, the 
intermediate grey represents the two-photon contribution and the light grey indicates the dominant four- 
fermion background) and the points with error bars are the OPAL data. The dashed empty histograms 
show a sample signal with (M, AM) =(100 GeV, 76.9 GeV), the empty histograms with the solid line 
show a sample signal with (M, AM) =(100 GeV, 49.2 GeV), and the dotted empty histograms show a 
sample signal with (M, AM) =(100 GeV, 12.2 GeV), all with arbitrary normalization. 
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Figure 2: Sample distributions of Lr for Standard Model background in the hadronic channel (solid 
histogram), chargino signal with arbitrary normalization (dashed histogram) and OPAL data at = 
207 GeV (points with error bars) at (M, AM) = (100 GeV, 49 GeV) and (100 GeV, 77 GeV). 

are less powerful discriminants or because they are too strongly correlated with better 
variables that are already in the likelihood. 

The following variables are used in the likelihood. Distributions of these variables are 
shown at preselection level in Figure 3. 



• M vis / \/s discriminates against four-fermion background, especially Wez/ 
with a fake lepton in one of the two jets; 



qqev 



• M'miss/v^ helps to reject W + W — > qqiu; 

• the scaled magnitude of the momentum of the lepton candidate, pt/yfs, also helps 
to reject W + W~ — > qqlv ; 

• M had /y/s also helps to reject W + W~ events; 

• or., the magnitude of the component of the missing momentum perpendicular to 
the event thrust axis, scaled to the beam energy, discriminates against two-photon 
and r + r~ events; 

• when the event, excluding the lepton candidate, is forced into two jets, a co P , the 
acoplanarity angle of the two jets, rejects two-photon backgrounds and events with 
ISR; 

• | cosset |, where 9j et is the polar angle of the jet closer to the beam axis, rejects 
two-photon and two-fermion backgrounds; 

• the scaled energy of the less energetic jet, Ej/y/s, helps to remove some two-fermion 
backgrounds and events with poorly reconstructed jets. 
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The following variables are used only for additional mass-dependent cuts. 

• the energy of the more energetic jet, scaled to the beam energy, discriminates against 
four-fermion backgrounds; 

• p™ lss / y/s, discriminates against two-photon backgrounds; 

• | cos^ a |, where 9 a is the polar angle of the a>T vector, rejects two-photon and t + t~ 
backgrounds; 

• | cos# miss | rejects two-photon and ISR events. 

Signal efficiencies for semileptonic events after these additional cuts, including the cuts 
on the likelihood variables, vary from about 0.08% at AM = 3 GeV to about 30% for 
5 GeV < AM < 10 GeV to about 60% for AM > 10 GeV. The main background 
surviving after the cuts in the region where AM > 10 GeV comes from four-fermion 
processes with a cross-section of about 1 - 2.5 pb. About 0.5 - 0.7 pb of e + e~T + T~ 
two-photon events survive for all values of AM and this is the dominant background for 
AM < 10 GeV. There are also about 0.05 - 0.1 pb of e + e~qq two-photon events for 
all AM and about 0.1 pb each of tau pairs and qq(7) events for AM > 10 GeV. Some 
examples of likelihood distributions are shown in Figure 4 for signal, data and expected 
Standard Model background. 

3.5 Neutralino Associated Production Hadronic Selection 

The final state topology for X2X1 production with a hadronic decay of the virtual Z° 
contains two hadronic jets and missing energy from the two stable lightest neutralinos 
in the final state. The analysis is very similar to the hadronic decay channel analysis 
for chargino pair production. Some changes are made to reflect the difference in event 
topology. The following variables are used: 

• Cut Variables — as for the chargino analysis except that no cuts are made on the 
number of jets in the event, on Egi nor on £j e ti; 

• Likelihood Variables: 

— |cos0 miss |, p™ lss , acO p and M vis are used, as for the hadronic channel of the 
chargino analysis; 

— -R™ d , the ratio of energy deposited in the forward region, | cos#| > 0.9, to the 
total visible energy, helps to remove ISR and two-photon events; 

— M vis /E vis , the ratio of the visible invariant mass to the total visible energy 
gives a measure of the "softness" and helps to remove qqz/z/ events; 

— <i 2 3 = Jy2?,El- ls) helps to remove backgrounds containing more than two jets. 
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Figure 3: Comparison of data and Monte Carlo distributions of selected likelihood variables for the 
semileptonic chargino analysis, shown at preselection level for ^/s = 207 GeV. The solid histograms 
represent the expected background, with dark grey for two-fermion events, cross-hatching for four-fermion 
events and single-hatching for two-photon events. The points with error bars are the OPAL data taken in 
the year 2000 at an average centre-of-mass energy of 206.1 GeV. Dotted histograms show a "W + W~-like" 
signal with (M, AM) =(80 GeV, 76.3 GeV) and dashed histograms show a "two-photon-like" signal with 
(M, AM) =(85 GeV, 6.9 GeV), both with arbitrary normalization. 
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Figure 4: Sample distributions of Lr for Standard Model background in the semileptonic channel (solid 
histogram), chargino signal with arbitrary normalization (dashed histogram) and OPAL data at y/s = 
207 GeV (points with error bars) at (M, AM) = (76 GeV, 70 GeV) and (100 GeV, 48 GeV). 

4 Extended Likelihood Method, Combination of Chan- 
nels and Energies 

Once the discriminant value Lr has been calculated for all events, as described in Sec- 
tion 3, it can be used either in an optimized cut analysis or in an "extended" likelihood 
analysis, as described in [8]. The optimized cut analysis has the disadvantage of being 
difficult to use effectively in a search such as this one, where the signal cross-section and 
hence the signal purity are unknown, and the number of selected events is sensitive to 
the position of the cut so that small modifications to the analysis can occasionally have 
a very large effect on the limits obtained. 

In the case of an unknown signal cross-section, an "extended" likelihood technique is 
more appropriate, as it is optimal for all luminosities and signal cross-sections. In this 
case, the distributions of the discriminant Lr are evaluated for signal and background 
samples and the value of Lr is evaluated for each data event. A single likelihood L can 
then be constructed which describes the probability, given the signal and background Lr 
distributions S and B, of observing N events passing all cuts with likelihoods Lr 4 : 

-u N N 

L-^rUPiL^S) (1) 

i=l 

where v = /i# + e£u;o",s is the expected number of background events passing the cuts, 
plus the expected number of preselected signal events for a signal cross-section as, and 
integrated luminosity C, assuming signal efficiency e. In the extended likelihood method, 
efficiencies are evaluated at the end of the cut-based selection. All candidates passing 
the preselection and additional cuts contribute to the final limit calculation, with weights 
given by their likelihood of being signal. There is no cut on the likelihood value, and 
therefore no sensitivity to cut position, making the limit very robust. 
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A weighting factor, u = 7 m is included to allow the cross-section limit to be set 
at a particular centre-of-mass energy using data taken at many different energies. Since 
charginos and neutralinos are fermions, a is assumed to scale as (5/s. With some algebraic 
manipulation [8] Equation 1 can be rewritten in the form: 

N 

\nL(a s ) = ~{ii B + eCuja s ) + Y J H^BB{L Ri )+eCuja s S{L Ri )} (2) 

i=i 

This can be calculated separately for each channel with data taken at several centre- 
of-mass energies. The 95% confidence level limit on the signal cross-section (Js (y/snmit) 
at ^/siimit = 208 GeV is computed using Bayes' theorem with a uniform prior in as by 
solving: 

/ >CT < 1 ) i ™ it v 

/ T-[ JV channels, ,/S r / —limits J —limit 

/ llj=l ^i{ a S ) aa S 

0-95 = J ° N (3) 

/ TT JV channels, VS r / limit N A -.limit 

/ lli=i L/iKPs ) aa S 

Jo 

Limits on the signal cross-sections are calculated according to Equation 3 separately 
for each channel in the chargino analysis and for all channels combined assuming 100% 
W* branching ratios. In all cases, all the energy bins from Table 2 are combined. Data 
from a particular centre-of-mass energy contribute to the combination only up to the 
kinematic cut-off for that energy, about (v^ - 1 GeV)/2. These cut-offs are then rounded 
down to the nearest 0.5 GeV. Thus, all data with energies above 192 GeV contribute up 
to a chargino mass of 95.5 GeV, only those with energies above 196 GeV contribute up to 
97.5 GeV and so on until the limit for chargino masses between 103 and 103.5 GeV comes 
only from the data taken at an average centre-of-mass energy of 208 GeV. The limit is 
calculated on the cross-section at ^/sumit — 208 GeV. 

For points where the background likelihood histograms contain empty bins, the ex- 
tended likelihood method cannot be used and a simple Poisson probability is evaluated, 
based on the observed number of events passing the preselection and additional cuts and 
the number expected from background Monte Carlo. This happens for most of the points 
with AM < 15 GeV and for some points with AM between 15 and 30 GeV. 



5 Systematic Errors 

There are numerous possible sources of systematic errors for these analyses. In general, 
however, it is found that even large errors on the efficiencies and backgrounds have rather 
a small effect on the limits, so only the larger effects are considered. Some sources of 
potentially large systematic uncertainties are: 

• limited Monte Carlo statistics for signal efficiency evaluation; 

• limited Monte Carlo statistics for background evaluation; 

• mismodelling of likelihood variables from incomplete or incorrect Monte Carlo sim- 
ulation and detector calibration effects; this includes, among other things, the effect 
of ignoring interference between multiperipheral and other four-fermion processes 
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containing an electron and a positron in the final state, for which a specific check 
was also done, and the effect was found to be negligible for the background, given 
the statistical errors on the multiperipheral background samples; a relative error of 
5 % was assigned to the signal efficiencies; 

• a further 1 % relative systematic error was assigned to the signal efficiencies to 
account for the uncertainty in the smoothing procedure applied to the reference 
histograms; 

• interpolation of signal histograms to obtain a finer grid in AM than was generated; 
a test was done using Monte Carlo samples generated at intermediate values of AM 
and a relative error of 1 % was assigned to the efficiencies; 

• effects of having signal generated only at three y/s to evaluate the efficiency at all 
energies and using rather crude (3/s models to evaluate the cross-sections of these 
fermionic signals at other ^fs were evaluated by repeating the whole analysis using 
reference histograms created from signals and backgrounds generated at 196 GeV 
instead of 207 GeV (in the kinematic range where this was possible) and comparing 
the results with the standard analysis; a conservative relative error of 5 % was 
assigned to the efficiencies; 

• uncertainties in the estimation of the integrated luminosity of about 0.2-0.3%. 

The evaluation of the uncertainties due to the first three items in the list is described in 
detail in subsections 5.1-5.2. Detector efficiency losses due to occupancy or due to the 
cuts on energy in the forward detectors, which reduce the effective integrated luminosity 
by about 3% were included simply by reducing all integrated luminosities by 3%. The 
following potential sources of systematic uncertainty were also considered and found to 
be negligible compared with those listed above: 

• effects of having two-photon background samples generated only at some of the \fs 
used and simply assuming that the cross-sections of these processes are proportional 
to log s to evaluate them at other energies. 

• effects of using reference histograms generated at a single ^/s for evaluating the like- 
lihood at all energies; a cross-check was done using reference histograms generated 
at = 196 GeV; 

All errors are evaluated by estimating the uncertainties on the rates of events passing 
a cut, since evaluating changes in the shapes of the Lr distributions at every point on the 
(M, AM) grid for every possible systematic is not a practical approach. 

5.1 Statistical uncertainties 

The statistical uncertainties on the signal efficiency and on the background estimates due 
to the finite number of Monte Carlo events generated are evaluated after the preselection 
and additional cuts for each mass grid point. 
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The most important sources of background are four-fermion final states, qq(7) and 
t + t~ two-fermion final states, and two-photon processes, principally e + e~r + r~ and un- 
tagged e + e~qq. The statistical errors on the two- and four-fermion final states are very- 
small, since the integrated luminosity of the generated Monte Carlo is about 50 to 100 
times the integrated data luminosity at most centre-of-mass energies; however, the two- 
photon processes have very large cross-sections and it was only possible to generate about 
8 to 15 times the integrated data luminosity. The statistical errors on each background 
are then added in quadrature to obtain the statistical error on the total background at 
each mass grid point for each centre-of-mass energy and in each analysis channel. The 
statistical errors on the background are therefore relatively high in the low AM region, 
which is dominated by two-photon processes, and almost negligible for AM > 10 GeV. 

5.2 Uncertainties due to variable mismodelling 

Systematic uncertainties due to differences between the data from the OPAL detector 
and the Monte Carlo models are estimated after an optimized cut on Lr. The Lr cut is 
optimized by selecting the bin boundary on the Lr distributions for signal and background 
which gives the smallest expected value of N$ 5 (the most signal events there could be, 
at a confidence level of 95%, given the expected signal and background distributions 
and assuming a number of observed events equal to the expected number of background 
events). The cut is optimized for the default Lr distributions, and the number of events 
passing it is evaluated with modified L R distributions made from reference histograms 
where variables have been shifted or smeared, to determine the systematic effects on the 
expected background levels if the background Monte Carlo differed from the data. 

For the evaluation of systematic effects due to variable mismodelling, reference his- 
tograms of the likelihood variables are constructed after the general preselection, with 
none of the additional mass-dependent cuts. Histograms for the sum of all the expected 
backgrounds are compared with histograms of the data. The peaks in the histograms 
are fitted with Gaussians for both the data and the Monte Carlo. The differences be- 
tween the peaks are compared with the statistical uncertainties on each background. If 
a distribution has two peaks and one is due mainly to two-photon events and the other 
mainly to four-fermion events, which can be determined by looking at the Monte Carlo 
data for each background individually, the comparison for each background is made only 
for the relevant peak. No attempt was made to apply this procedure to the variables with 
obviously non-Gaussian distributions, such as cosine distributions; their shapes generally 
show good agreement between data and Monte Carlo. If the difference in peak position 
is larger than the statistical uncertainty for a particular background, then the difference 
minus the statistical uncertainty is taken to be the systematic uncertainty, and reference 
histograms are constructed for each background with the variable shifted by the amount 
given by this prescription. A similar procedure is applied to the comparison of the widths 
of the peaks in data and Monte Carlo and used to determine an amount by which to 
smear the variables for the reference histograms. When these are compared with quanti- 
ties obtained by comparing OPAL Z° data with Monte Carlo, which allows a comparison 
which is not statistically limited, agreement is reasonably good. 
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Table 4: Relative systematic errors on major backgrounds in the scmilcptonic channel for the chargino 
analysis. 



Background 


Relative 


Source 


Systematic Error 


qq 


0.48 


r + T - 


0.19 


4-fermion 


0.25 


e + e~ 


1.0 


e + e~r + T~ 


0.95 


e + e~qq 


1.0 



Table 5: Relative systematic errors on major backgrounds in the hadronic channel for the chargino and 
neutralino analyses. 



Background 


Relative 


Source 


Systematic Error 


qq 


0.3 


4-fermion 


0.8 


e + e~qq 


0.5 



5.3 Method for extrapolating systematic uncertainty estimates 
to the full mass grid 

Variations due to systematic uncertainties are evaluated at a centre-of-mass energy of 
207 GeV, which allows mass grid points near the kinematic limit to be considered, and 
is an energy at which both signal and background Monte Carlo samples are available. 
It is assumed that the uncertainties at 207 GeV are typical and can be used at all 
centre-of-mass energies. The errors are evaluated at three mass grid points: (M x ± = 
80.0 GeV, AM = 76.3 GeV), where signal events give a final state topology almost identi- 
cal to that of on-shell W-pair production, (M x ± = 85.0 GeV, AM = 6.9 GeV), where the 
detected particles from signal events are very soft and the main background is from two- 
photon events, and (M x ± = 100.0 GeV, AM = 49.2 GeV), which is near the kinematic 
limit and has background coming from both four-fermion and two-fermion production. 
Uncertainties are estimated for each of the main background sources at each of these 
points. The difference between the result with the smear or shift and the original result 
is found and the statistical uncertainty for that background subtracted off. These differ- 
ences are then summed in quadrature for the shifts and smears of all the variables. The 
errors at the three points are then compared and, in general, the error from the most 
relevant of the three for the particular background is chosen to be the error; in other 
words, the errors on the four-fermion background are generally taken from the point 
(M x ± = 80.0 GeV, AM = 76.3 GeV), while the errors on two-photon backgrounds come 
from the point (M x ± = 85.0 GeV, AM = 6.9 GeV). The results of this study are used to 
give the relative systematic errors shown in Tables 4 and 5. The fractional contribution of 
each background was then calculated for the rest of the mass grid and the relative errors 
on each background scaled according to the contribution of that background at each point. 
The systematic errors for each background at each grid point are then added in quadra- 
ture and eventually summed in quadrature with the statistical error on the background. 
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The efficiency errors are summed in quadrature and a total relative systematic error of 
7% is applied to the signal efficiencies and summed in quadrature with the statistical 
uncertainties. 



5.4 Incorporating uncertainties into cross-section limits 

Once the systematic errors have been evaluated for the signal efficiency and background 
expectation at every point on the mass grid for all energies, the limit setting code is 
run twice more over the data: once with all the efficiencies decreased by their errors 
and all background expectations increased by their errors to give the worst expected 
sensitivity, and once with the efficiencies increased and backgrounds decreased to give 
the best expected sensitivity. There are then three limits for each point in the grid: the 
standard one, L95q, the one standard deviation "worst case" one, L95+, and the one 
standard deviation "best case" one, L95_. All are evaluated at a centre-of-mass energy 
of 208 GeV. A Gaussian-weighted average of these gives the limit convolved with the 
systematic errors: 

L95 = ^ ^ -. (4) 

1 + 4= 

This is a three-point integral; it is a simple and robust way to deal with large errors, 
which tend to make more complex methods unstable. Despite the large relative errors 
on the background expectation, the limits typically change by only about 0.001-0.02 pb, 
with an increase of about 0.1 pb for the hadronic channel in the "W + W~" region where 
(M, AM) « (80 GeV, 80 GeV) and an increase of several picobarns along the AM = 
3 GeV "two-photon" strip for the semileptonic channel. 



6 Final cross-section limits 



In order to give a visual impression of what the limits might look like in the absence of 
a signal, "expected" limits are calculated. For each point on the mass grid, an average 
is made of the cross-section limits from a hundred toy Monte Carlo experiments. In 
each experiment, a random number of "events" is generated using a Poisson distribution 
with a mean value of hb, the expected number of background events for the mass grid 
point. Each "event" is assigned a random likelihood from the background likelihood 
distribution B(L Ri ). The signal cross-section limit is then calculated just as it is for the 
data. 



6.1 Chargino Pair Production, Hadronic channel 

The limit observed is in general slightly higher than the expected limit; however, the 
95% confidence level limit on a^+^-BR 2 (xf — > qqx?) (shown in Figure 5) is below about 
0.3 pb almost up to the kinematic limit for AM < 3 GeV and well below 0.1 pb in most 
of the parameter space explored, whereas a chargino signal would be expected to have a 
cross-section of several picobarns. 
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Expected limit at 208 GeV, Hadronic Channel Observed limit at 208 GeV, Hadronic Channel 




M(xD (GeV) MOcD (GeV) 

Figure 5: The expected and observed 95% confidence level limits on a(e + e~ — > XiXi^^ix^ ~> QQXi)i 
using data taken at 192-209 GeV. White areas indicate that the limit is greater than 2 pb or that no 
limit could be set. 

6.2 Chargino Pair Production, Semileptonic channel 

The observed limits are not as stringent as the expected ones; however, the 95% confidence 
level limit on a^±^±BR{xi q.QXi)BR(xf ^ u Xi) (shown in Figure 6) is below about 
0.3 pb almost up to the kinematic limit for AM < 3 GeV and well below 0.1 pb in much 
of the parameter space explored. 



Expected limit at 208 GeV, Semileptonic Channel Observed limit at 208 GeV, Semileptonic Channel 




MOcD (GeV) MOcD (GeV) 

Figure 6: The expected and observed 95% confidence level limits on the cross-section for the final state 
containing jets and a lepton with missing energy, using data taken at 192-209 GeV. White areas indicate 
that the limit is greater than 2 pb or that no limit could be set. 



6.3 Chargino Pair Production, Leptonic channel 

The observed 95% confidence level limits on er~± ^±BR 2 (xt ~^ ^ u Xi) taken from [4] are 
shown in Figure 7 and are slightly better than the expected limits. 
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Expected limit at 208 GeV, Leptonic Channel Observed limit at 208 GeV, Leptonic Channel 




MOcD (GeV) MOcD (GeV) 

Figure 7: The expected and observed 95% confidence level limits on the product of the cross-section for 
chargino pair production and the square of the branching ratio for x^ 1 to leptons and missing energy, 
using data taken at 192-209 GeV. White areas indicate that the limit is greater than 2 pb or that no 
limit could be set. 

6.4 Chargino Pair Production, Channels Combined Assuming 
100% Branching Ratios 

If the lightest neutralino is the LSP, chargino branching ratios to final states containing 
qq and tv are the W ± decay branching ratios, except for very low values of AM where the 
leptonic branching ratio increases rapidly. The three channels are combined according to 
W branching ratios and these results, which are generally valid in the absence of light 
scalar leptons and scalar neutrinos, are shown in Figure 8. 



Expec ted limit at 208 GeV, All Channels, W * BR Obser ved limit at 208 GeV, All Channels, W * BR 




M(xD (GeV) M(xD (GeV) 



Figure 8: The expected and observed 95% confidence level limits on the cross-section for chargino pair 
production, assuming 100% W* branching ratios for the chargino decays and using all the data sets 
analyzed. White areas indicate that the limit is greater than 2 pb or that no limit could be set. 
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6.5 Neutralino Associated Production, Hadronic Channel 



Since only the hadronic channel was analyzed for neutralino associated production, the 
results are presented both as a limit on a^oBRijQ qqx?) (Figure 9) and as a limit on 
(j^o^,o assuming 100% Z°* branching ratios for x° decay (Figure 10). The 95% confidence 
level limit on the cross- sect ion is less than 0.1 pb for AM > 5 GeV to within about a 
GeV of the kinematic limit. There is no value of (M^), M(xi)) for which a significant 
excess is observed. 



Expected limit at 208 GeV, Hadronic Channel Observed limit at 208 GeV, Hadronic Channel 




Figure 9: The expected and observed 95% confidence level limits on the product of the cross-section for 
X2X1 associated production and the branching ratio for - * lOXi decays. No signal Monte Carlo events 
were generated with M^o + M^o < 100 GeV, so no limit is evaluated in this region. White areas indicate 
that the limit is greater than 0.5 pb or that no limit could be set. 



Expected limit at 208 GeV, Hadronic Channel, Z° BR Observed limit at 208 GeV, Hadronic Channel, Z° BR 




Figure 10: The expected and observed 95% confidence level limits on the cross-section for X2X1 associated 
production, based on decays to hadronic final states and assuming 100% Z°* branching ratios for the X2 
decays. No signal Monte Carlo events were generated with M^o + M^o < 100 GeV, so no limit is evaluated 
in this region. White areas indicate that the limit is greater than 0.5 pb or that no limit could be set. 
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Despite the small general excess which results in larger than expected cross-section 
limits in the hadronic and semileptonic channels, there is no indication of chargino or 
neutralino production in the OPAL data. 

7 Interpretation in the Constrained MSSM 

The experimental information required to obtain limits in any particular model is con- 
tained in the cross-section limit for chargino production as a function of the leptonic 
branching ratio and in the cross-section limit for neutralino associated production as a 
function of the invisible branching ratio. Comparing the limits with the cross-section and 
leptonic branching fraction predicted by any specific model allows that model to be tested 
with the OPAL data. 

In this section, the results are interpreted in the context of the Contrained Minimal 
Super symmetric Standard Model (CMSSM). This is a model inspired by Minimal Super- 
gravity, in which it is assumed that at some large Grand Unification (GUT) scale all of the 
gauginos have a common mass m^. This implies that there is a fixed relation between 
the three gaugino masses at the electroweak scale, M\ for U(l) partners, M 2 for SU(2) 
partners, and M 3 for SU (3) partners, so a single free parameter, taken to be M 2 , generates 
all the gaugino masses. It is also assumed that the sfermions (but not all scalars as in 
true Minimal Supergravity) have a common mass at the SUSY-breaking scale, denoted 
by itlq. With these constraints, there are only six free parameters not present in the SM: 
M 2 , mo, the common trilinear sfermion-Higgs coupling at the Planck scale A , the ratio 
of the two vacuum expectation values for the Higgs doublet fields coupling to down- and 
up-type quarks tan (3, the mass mixing parameter of the Higgs fields /z, and the mass of 
the pseudoscalar Higgs boson ttia- With these six parameters, it is possible to calculate 
cross-sections for X1X1 an d X2X1 production. 

The properties of the charginos and neutralinos depend chiefly on M 2 , /i and tan/5. 
If M 2 3> then the lightest neutralino is higgsino-like and t-channel contributions are 
negligible; the xt , Xi an d X2 are au close in mass and the cross-section for X2X1 production 
is similar to that of X1X1 production. If M 2 <C \fJ>\ then the Xi is gaugino-like; t-channel 
processes are important and may interfere destructively with the s-channel. Scans over 
M 2 and /i are done at several values of tan/5. The m and A parameters determine 
the masses of the sfermions. Results here are shown for m = 500 GeV, implying that 
there are no light sleptons or sneutrinos, and two-body decays of charginos to £ ± z> can 
be ignored, and for A = 0, which generally means that mixing of left- and right-handed 
sfermion eigenstates is assumed to be small. A scan is performed with SUSYGEN [33] over 
the region < M 2 < 5000 GeV, -1000 GeV < fi< 1000 GeV at tan/3 = 1,2,5, 10,40. 
The granularity of the scan depends on the region and is generally between 1 and 5 GeV 
in M 2 and /i. 

The simplest thing to do with the cross-section limits obtained in the previous section 
and the results of the CMSSM scan is to see what regions of the parameter space are 
excluded. By construction, the three chargino channels do not use the same events, so 
the best way to use this information is to sum the limits for the three channels according 
to the branching ratios predicted by SUSYGEN for a particular set of M 2 , //, tan (3 values 
and compare with the SUSYGEN prediction for the chargino total cross-section. The 
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events selected in the neutralino associated production search may overlap with the events 
selected in the chargino search, so the optimal way to include both searches is to calculate 
whether the expected limit for the chargino or the neutralino search would give the better 
exclusion for a particular M 2 , fi, tan /3, and use the channel which is expected to be better. 
Exclusions in the M 2 ,/i plane are shown in figure 11 for the five values of tan/3 studied. 
These exclusions are valid only for m > 500 GeV and for A = 0. The regions which 
are shown as kinematically accessible but not excluded correspond to regions where X2X1 
production is possible, but with an expected cross-section too small to be excluded. 

The final step is to see what masses of charginos and neutralinos are excluded in the 
CMSSM. The kinematically accessible regions of the M-±,M^o and M^o,M^o planes are 
divided into a 1 by 1 GeV grid. All of the scanned (M 2 ,/x,tan/3) points are examined 
to see what masses they predict for the lightest chargino and two lightest neutralinos. 
The smallest cross-section predicted in any 1 by 1 GeV grid square is compared with the 
cross-section limit. The CMSSM exclusions for the chargino mass plane are shown for 
the five values of tan/3 in figure 12 and the neutralino mass plane exclusions are shown 
in figure 13. 

Assuming that m > 500 GeV and A = 0, the following conclusions can be made. 
For AM > 5 GeV, charginos with masses less than 101 GeV are excluded for all tan/3. 
For AM > 5 GeV, x° with masses less than 78 GeV and x? with masses less than 40 GeV 
are excluded for all tan (3. 

8 Conclusions 

No hint of charginos or neutralinos was observed in the data from OPAL. The cross- 
section limits presented in this paper for charginos or neutralinos decaying to particular 
final states may be regarded as completely model-independent and can be used to set 
limits on any model. Limits are set on the cross-section for chargino pair production by 
combining the three chargino final states under the assumption that all charginos decay 
to a virtual W, which then decays with the usual branching ratios to leptons and hadrons. 
Limits are calculated for the specific example of the Constrained MSSM in the absence 
of light sfermions or mixing between right- and left-handed sfermions. Charginos are 
excluded almost up to the kinematic limit set by the maximum LEP beam energy of 
about 104 GeV. 
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Figure 11: Regions of the CMSSM parameter space excluded by the OPAL search are shown in black. 
These include regions excluded by the OPAL data taken at energies up to 189 GeV [5] as well as those 
excluded with the present dataset and analysis. The regions with the lighter shading are kinematically 
accessible for either xtxl or X2X1 production at y/s = 208 GeV, but are not excluded. These results 
are valid for mo > 500 GeV and Ao = and for the values of tan /3 indicated on the individual plots. 
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Figure 12: The shaded regions of the xt — Xi mass plane are excluded by OPAL data, with the lighter 
colour showing the region excluded by previously published searches [5] at centre-of-mass energies up to 
189 GeV, and the black region being excluded by the present search. The kinematic limits of the chargino 
search are indicated by dotted lines; the very small excluded regions outside these limits are due to the 
neutralino search. The cross-hatched regions labelled "No CMSSM solution" correspond to regions of the 
parameter space having no solution in the CMSSM. The regions left white are not excluded. 
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Figure 13: The shaded regions of the %° — Xi mass plane are excluded by OPAL data, with the lighter 
colour showing the region excluded by previously published searches [5] at centre-of-mass energies up 
to 189 GeV, and the black region being excluded by the present search. The kinematic limits of the 
neutralino search are indicated by dashed lines. The dashed line at M^o + M%a = 100 GeV is simply a 
reminder that for the present analysis, no X2X1 Monte Carlo events were generated with lower masses; 
exclusion shown in black below this line is purely due to the chargino searches. The singly hatched region 
labelled "LEP 1 Search" is the region kinematically accessible to LEP 1 searches; it was not considered 
in this study. The cross-hatched region labelled "No CMSSM solution" corresponds to regions of the 
parameter space having no solution in the CMSSM. The regions left white cannot be excluded by the 
data. 
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